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ABSTRACT 

It is commonly believed that the earliest stages of star-formation in the Universe were 
self-regulated by global radiation backgrounds - either by the ultraviolet Lyman- 
Werner (LW) photons emitted by the first stars (directly photodissociating H 2 ), or by 
the X-rays produced by accretion onto the black hole (BH) remnants of these stars 
(heating the gas but catalyzing H2 formation). Recent studies have suggested that a 
significant fraction of the first stars may have had low masses (a few M Q ). Such stars 
do not leave BH remnants and they have softer spectra, with copious infrared (IR) 
radiation at photon energies ~ leV. Similar to LW and X-ray photons, these photons 
have a mean-free path comparable to the Hubble distance, building up an early IR 
background. Here we show that if soft-spectrum stars, with masses of a few M Q , 
contributed > 1% of the UV background (or their mass fraction exceeded ~ 90%), then 
their IR radiation dominated radiative feedback in the early Universe. The feedback is 
different from the UV feedback from high-mass stars, and occurs through the photo- 
detachment of H~ ions, necessary for efficient H2 formation. Nevertheless, we find 
that the baryon fraction which must be incorporated into low-mass stars in order to 
suppress H2-cooling is only a factor of few higher than for high-mass stars. 

Key words: cosmology: theory - early Universe - galaxies: formation - molecular 
processes - stars: Population III 



• . 1 INTRODUCTION 

3 . 

^ In hierarchical models of structure formation, the first 
stars in the Universe form in dark matter (DM) 
minihalos with masses of ~ 10 5 M© at redshifts of 
z ~ 20 — 30, through eff i cient cooling of the gas by 
H2 (iHaiman. Thoul fe Loebl 1 19961 : see a comprehensive 
review by iBarkana fc Loebl 12001! ). However, soon after 
the first stars appear, early radiation backgrounds begin 
to build up, resulting in feedback on star-formation. In 
particular, the UV radiation in the Lyman-Werner (LW) 
bands of H2 can photodissociate these molecules and 
suppress gas cooling, p ossibly p reventing star-fo r mation 
jHaiman. Rees fc Loebl 1 19971: lOmukai fc Nishil 1 1999 1; 



the early minihalos were very massive (~ 100 Mq), ow- 
ing to the rapid mass accretion enabled by H2 cooling. 
These stars wo uld then leave b ehind remnant BHs with 
similar masses (|Heeer et al.ll2003h . and produce X-rays, ei- 
ther by direct accretion or by forming high-mass X-ray 
binaries. A soft X-ray background at photon energies of 
> IkeV, at which the early intergalactic medium (IGM) 
is optically thin, then provides further global feedback: 
both by heating the IG M, and by catalyzing H2 fo r ma- 
tion i n collapsing halos ijHaiman. Rees fc Loebl 19961: 
200 it IVenkatesan, Giroux fc Shulll l200ll: iGlover fc Brand! 



20031; iMadau et al.l |2004 IChen fc Miralda-Escudel 12004 
Ricotti. Ostriker fc Gned in 2005; Mirabel et al.ll201ll 



Haiman. Abel fc Rees! boOCt ICiardi. Ferrara fc Abell 1200" 
Machacek Bryan fc Abell l200ll: iRicotti. Gnedin fc Shulj 
200l| . 120021: iMesinger. Bryan fc Haimanll2006l; IWise fc Abel 
20071 : lO'Shea fc Normanl 120081 : Ijohnson. Greif fc Brornin 
Abell l2008allb| : IWhalen et all 120081 : 



Wise 



Mesi nger. Bryan fc Haim an 2009) 



Numerical simulations (e . g. Abel. Bryan fc Normanl 



120021 ; iBromm. Coppi fc Larsonl |2002| : lYoshida et al.l 120031 ) 
have long suggested that the metal-free stars forming in 
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Recent simulations have been pushed to higher spa- 
tial resolution, and in some cases, using sink particles, 
were able to continue their runs beyond the point at 
which the first ultra-dense clump developed. The gas in 
the central regions of at least some of the early minihalos 
were found to fragment into two or more distinct clumps 
llTurk. Abel fc O'Sheal 120091: IStacy. Greif fc Brornml l2010l : 



iGreif et al.ll201ll : IClark et alll201ll : iPrieto et alj|201lf ). This 
raises the possibility that the first stars formed in multiple 
systems, and that many of t hese stars had low er masses than 
previously thought (but see iTurk et alj|2012l for still higher 
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resolution simulations that suggest less efficient fragmenta- 
tion). 

There is also some observational evidence suggesting 
a lower characteristic Pop III mass. Massive (> 140 Mg), 
non-rotating metal-free stars are expected to end their 
lives as pair- instability supernovae (PISNe), and the non- 
detection of the characteristic PISN nucleosynthetic pat- 
terns in metal-poor stars suggests that the typical Pop 
III stars did not form with such high masses (see, e.g., a 
recent review by iFrebel fc Norrisl 1201 lh . The observations 
of carbon-enhanced metal poor stars may further imply a 
significant numbe r of Pop HI stars with masses as low as 
M = 1 - 8 M (|Tumlinsonll2007bl lal). Finally, the recent 
discoveries of extremely metal poor stars with no sign of 
C or N enhancement shows that low-mass star formation 
could o ccur at metallicitie s much lower than previously as- 
sumed llCaffau et alj|201ll. 12012:), likely facilitated by dust 
fragmentation ( Schneider et al.ll2012h . 

Motivated by the above, in this Letter, we examine ra- 
diative feedback from an early cosmic IR background, pro- 
duced by a population of low-mass stars. Although we fo- 
cus on low-mass PopIII (i.e., metal-free) stars, our conclu- 
sions are more generic, and apply at any cosmic epoch when 
significant numbers of low-mass PopII stars co-exist with 
massive PopIII stars. Low-mass stars a re expected to have 
soft spectra, even in the metal-free case ()Tumlinson fc Shulll 
l2000l : IMarigo et al.ll200ll ; ISchaererl [20021 ). producing signifi- 
cant radiation at ~ leV, near the photo-detachment thresh- 
old (0.76 eV) of the H~ ion. H~ is a reactant in the dom- 
inant formation channel for H2, (H~ + H — > H2 + e~) and 
its destruction can therefore have a dramatic impact on the 
thermal evolution of metal-free gas. 

In fact, it is well known that photo-detachment of H~ 
by cosmic microwave background (CMB) photons kept the 
H2 formation rate in the early Universe very low, until 
the CMB photons redshifted to l o wer energies at z ~ 100 
(Hira sawa. Aizu fc Taketanil Tl969: Hir ata fc Padmanabhanl 
I2006T I. H~ photo-detachment can become globally important 
again once stars begin to form, if they have soft spectr4E 

The aim of this Letter is to quantify (i) if and when, due 
to low-mass stars, H~ photo-detachment again became the 
dominant process to limit H2 cooling in the earliest proto- 
galaxies, and (ii) to what extent this may have increased or 
decreased the net global negative radiative feedback in the 
early Universe. We focus on the importance of this negative 
feedback for minihalos (as opposed to the more massive ha- 
los that cool even in the absence of molecular hydrogen). 
In order to accomplish this, we perform "one-zone" calcula- 
tions, following the coupled chemical and thermal evolution 
of the gas in the presence of a cosmological radiation back- 
ground, including H~ photo-detachment by IR photons and 
H2-photodissociation by LW photons. 

The rest of this paper is organized as follows. In § [2] 
we describe our chemical and thermal modeling. In § [3] we 
present our results for the relative importance of IR radia- 
tive feedback, with various assumptions about the stellar 
populations; we also compare our results to previous stud- 
ies. Finally, in § [4] we offer our conclusions. Throughout 



1 To our knowledge, this po i nt wa s first noted and discussed by 
IChuzhov, Kuhlen fc Shapiro! J2007h : see § 13.31 below. 



this paper we adopt a standard ACDM cosmological back- 
ground model: (Rdm, »b, »a, h)= (0.233, 0.046, 0.721, 
0.701) (|Komatsu et alj|201ll ). 



2 MODELING 

2.1 Background Spectrum 

We assume that the early Universe is filled with back- 
ground radiation produced by stars. Massive metal-free 
stars have hard spectra, with effective blackbody temper- 
atures of T c ff ~ 10 5 K on the zero-age main sequence 
(ZAMS), nearly independent o f their mass above > 
100 Mm llMarigo et ail l200ll : iBromm. Kudritzki fc Loebi 
l200ll : ISchaereril2002l ). Below this mass, the effective temper- 
ature decreas es monotonically, and reaches T e g ~ 10 4 K at 
M+ « 2 M fl l|Tumlinson fc Shulll I2OO0I ; IMarigo et ai1l200ll : 
ISchaerenl2002l ). 

We first consider a composite spectrum produced by co- 
existing low- and high-mass stars, the relative abundances 
of which are allowed to vary. We choose characteristic masses 
of Mio = 1.2 M and M hi = 100 M , which have time- 
averaged effective temperature s of r e ff.io = 10 3 9 5 K and 
T eff ,hi = 10 4 88 K, respectively (jMarigo et al.ll200ll ). These 
stars have main sequence lifetimes of t ms ss 3 Gyr and 
~ 3 Myr. In our analysis below, for low-mass stars with 
tms < 0.5Gyr (the age of the universe at z » 10), we reduce 
the total radiative output by the factor (t ma /0.5Gyr). We pa- 
rameterize the population by the mass fraction of low-mass 
stars, /™ ass and also by their corresponding fractional con- 
tribution /i" v to the total UV radiation output (at 13.6eV). 

We caution that this is of course only a crude charac- 
terization of the true background - in principle, one needs 
to consider the time-evolving spectra of stars with a range 
of initial masses. Departures from a black-body shape may 
also be important, as our results are sensitive to photons 
with the particular energies of ~ 2eV (for H~ detach- 
ment) and ~ 12eV (for H2 dissociation). However, given 
that the stellar IMF is unknown, we here adopt this sim- 
ple prescription and defer a more detailed treatment to fu- 
ture work. In order to understand the dependence on stellar 
mass, in § 13. 21 we repeat our calculations assuming the back- 
ground is produced by stars with a single mass in the range 
0.7 M c h al ./M0 100 (corresponding to effective temper- 
atures in the range 6000 < T efi < 10 5 K). 

The early IGM is optically thin at the IR energies of 
~ 2eV, relevant for photo-detachment of H~ (owing to 
the very low abundance of both intergalactic H~ and of 
e~). However, before reionization, the IGM is opaque above 
13.6eV, and we accordingly assume zero flux above this en- 
ergy. UV photons in the LW bands (11.2-13.6eV), travel- 
ing over cosmological distances, will also be absorbed by HI 
once they redshift into resonance wi th a Lyman line. This 
resul ts in a "sawtooth modulation" l|Haiman. Rees fc Loebl 
1997), which reduces the H2 dissociation rate by about 
an order of m agnitude at z = 15 compared to the opti- 
cally thin rate (|Wolcott-Green fc Haimanfeoill ') . We include 
this reduction in our calculation. In principle, the optical 
depth in the LW lines due to intergalactic H2 itself, with 
an abundance of nm/wH ~ 10~ 6 , can be of order a few 
|Haiman. Abel fc Reedl200cl ; iRicotti. Gnedin fc Shulj|200ll : 
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iKuhlen fe Madaul 120051 ): this additional opacity could fur- 
ther reduce the ^-dissociation rate and would strengthen 
our conclusions; we conservatively ignore it in our calcula- 
tions. 



2.2 Chemistry and Cooling 

We model a static gas cloud which has condensed to 
the maximum density achievable by adiabatically collaps- 
ing in a minihalo with virial temperature T v i r : p ma x = 
/pPiGM(Tvir/TiGM) 3/2 - Here pi G M and Tigm are the den- 
sity and temperature of the smooth background IGM, and 
f p (chosen below) is of order unity and is used for calibra- 
tion against simulations. At this stage, only if sufficient H2 
forms will the cloud be able to radiatively cool and continue 
collapsing to higher densities, and ultimately form stars. 

We follow the chemical and thermal evolution for pri- 
mordial gas, using a standard chemical reaction network 
among nine species: H, H+, H~, He, He+, He 2+ , H 2 , H+ 
and e~ (and photo n s). R adiative cooling by H2 is modeled 



Galli fc Pallal (I1998T). Ou r chem ical model is adopted 



Shang, Bryan fc Haimanl (|2010h . where the reader is 



as m 
from 

referred for more details and references. We adopt initial 
conditions appropriate for a minihalo (i.e. with molecular 
hydrogen and electrons fractions of ih 2 — 2 x 10 -6 and 
x c — 10 -4 , and initial gas temperature T gas = T v i r = 400K). 
We follow the coupled thermal and chemical evolution using 
the stiff equation solver LSODAR for ~ 10 s years (or ~ 20 per 
cent of the Hubble time at z ~ 10, after which it is likely 
the halo will have merged with another to form a larger 
collapsed object). 

Photo-detachment of H~ by continuum photons 
(H~ + 7 — > H + e~) is dominated by ~ 2eV photons 
(at which the H~ detachment rate peaks for T e g w 
10 4 K, slightly above the 0.755eV thresh old). We fit the 
freq ue ncy-dependent cross-section as in IShapiro fc Kan j 
1 19871 ). convolving this with the blackbody spectrum 
to find the rate coefficient, parameterized as £25 = 
10~ 10 aJ 2 i s _1 . Here and below, J21 denotes the spe- 
cific intensity at the Lyman limit (13.6eV) in units of 
10~ 21 erg s -1 cm -2 sr _1 Hz -1 . For the composite spectrum 
with (T cff .i ,T cff , hi ) = (10 3 ' 95 K, 10 4 ' 88 K), we have a = 
aio/io V + Qhi(l - /D, and (a lo ,a hi ) = (10 4 ,0.17). 

For the photodissociation of H2 by LW photons we 
use the fitting formulae for the optically-thick rate pro- 
vided bv lWolcott-Green. Haiman fc Brvanl |201ll ). which in- 
cludes self-shielding as well as shielding by HI. The col- 
umn densities are specified by assuming the size of the col- 
lapsing region equals the Jeans length; however, we reduce 
-/Vh 2 by a factor of 10, which produces better agreement 
with the shielding found in three-dime nsional simulations 
l|Wolcott-Green. Haiman fc Brvaifeoill '). 



3 RESULTS AND DISCUSSION 
3.1 Composite Spectrum 

The Critical Flux 

To assess whether H2-cooling can be suppressed by the IR 
background, we first run one-zone models at various dif- 
ferent radiation intensities, normalized at the Lyman-limit 



10- 



1 .0 

s 



no H photo-detachment^- 
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Figure 1. Top panel: The critical flux J 2 r \ t required to suppress 
H2-cooling as a function of /j" v - Here is a proxy for the 
fraction of the Lyman limit flux contributed by low— mass stars 
with an effective temperature T c g ss 10 4 K (rather than high-mass 
stars with T e g rj 10 s K). The corresponding mass fraction in the 
low— mass stars, /j™ ass , is shown along the upper horizontal axis. 
Bottom panel: The fraction of baryons that must be converted 
into stars in order to produce the critical flux shown in the top 
panel (solid line). Also shown is the required fraction if H — photo- 
detachment is artificially switched off in the chemistry network 
(dashed line). 



J21 . We employ a Newton- Raphson scheme to find the "crit- 
ical flux," Jcrit, defined by requiring that the cooling time 
always remains longer than the dynamical time (eq. 10 in 
lHaiman. Rees fc Loebl[l997l ). Repeating this calculation for 
halos with different masses and redshifts, we have found 
that with the density normalization f p = 0.5, our resulting 
Jcrit = Jcrit (Afhaio, z) agrees well (to within a factor of two) 
with the values, and the m ass- and redshift-dependence of 
Jcrit found in simulations (Machacek. Bryan fc Abell 120011 ; 



iMesinger. Bryan fc Haiman 2006). Throughout this section, 
we show our results for a single halo mass scale, T v i r = 400K 
and redshift z = 18 (corresponding to the lowest-mass mini- 
halos that can cool via H2). Apart from a monotonic overall 
increase of J cr it with T v i r and with z, our results are qualita- 
tively the same for other redshift and halo masses (up to the 
scale corresponding to T v i r ~ 10 4 K, at which atomic cooling 
becomes important). 

Our main result is shown in Figure [1] the top panel of 
which shows the critical flux J cr it as a function of /i£ v (the 
corresponding mass fraction /!™ ass is shown along the upper 
horizontal axis). The critical flux decreases dramatically as 
/io V increases above the per cent level. Since the normal- 
ization of the critical flux is quoted at 13.6eV, close to the 
LW band, the critical LW dissociation rate would be nearly 
independent of f^; the decrease is caused entirely by the 
H~ detachment by IR photons. The IR radiation becomes 
more important than the LW radiation when the decrease 
reaches a factor of two, which occurs for /£ v » 10~ 2 . How- 
ever, because the low-mass stars emit far fewer Lyman-limit 
photons than those with effective temperatures ~ 10 5 K, a 



mass fraction of f" c 



0.9 is required to achieve even 
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this percent level contribution to the flux (as shown by the 
upper horizontal axis). For reference, the mass fraction of 
stars in the O.I-I.2M0 (O.I-2M0) range for the commonly 
used Chabrier IMF is 86% (94%). In the limit of purely low- 
mass _staxs 2 _ J 7crit is_decreased by a factor of w 150. (Note 
that lClark et al.ll201ll do find the protostellar IMF peaked 
at < IMq, but ZAMS masses are unknown.) 
Global Impact 

We next ask whether the critical flux J cr it=fewx (10~ 4 — 
xl0~ 2 ) can be produced by the mixture of high- and low- 
mass stars. To answer this, we determine the fraction of 
baryons that must be incorporated into stars to produce 
a given J cr it, as a function of /j" v . The energy density in 
LW photons u 7 ~ 47tz/ Jcrit /c is converted to a stellar mass 
density p* by computing the number of L W photons emitted 
per stellar baryon0 We use the data from lMarigo et al"1l200ll 
for this purpose, and find the mass-weighted average for the 
high- and low-mass populations. 

The bottom panel of Figure [T] shows the critical baryon 
fraction, /*, C rit, which must be incorporated into stars to 
achieve Jcrit (/i" v ) (normalized to its value in the absence of 
any low-mass stars). The fraction varies relatively little over 
the range ^ /ic, v ^ 1 (solid curve). The factor of ~ 10 in- 
crease as fiJ — ¥ 1 is due almost entirely to the long lifetimes 
of the low-mass stars, and the corresponding reduction of 
their radiation output over the finite 0.5 Gyr age of the uni- 
verse (at z ~ 10). The otherwise near- flatness of this curve 
is a coincidence: increasing decreases both the critical 
flux and the production of LW photons per stellar mass, and 
these two factors turn out to nearly cancel. The bottom-line 
is that H~ photo-detachment can indeed suppress subse- 
quent star formation; this requires converting baryons into 
low-mass stars up to an order of magnitude more efficiently 
than in the high-mass case. Finally, the dashed curve in 
the bottom panel in Figure [T] shows that if the IR photons 
from the low-mass stars were neglected, the required stel- 
lar density would increase by two orders of magnitude. This 
highlights the importance of including IR feedback in future 
models if the Pop III IMF indeed extends to such low masses 
(and also at all cosmic epochs when PopII stars are already 
present). 



3.2 Dependence on effective temperature 

How low do the masses of low-mass stars need to be before 
IR feedback becomes important? To answer this question, 
we next consider stellar populations with a single effective 
temperature, in the range 6000 < T e g/K < 10 5 , and com- 
pute how Jcrit and /», cr it depend on T e g. The top panel of 
Figure [2] shows a dramatic drop in the critical flux owing to 
the copious ~ 2eV photon output below T c g ~ 1.5 x 10 4 K, 
or M < 2M0 . As shown in the bottom panel of this figure, 
the density in stars required to produce J cr it remains essen- 
tially constant down to this T e g, and very near the value 
required for LW feedback by high-mass stars. This flatness 
again results from the cancellation of two effects: the criti- 
cal flux and the LW photon output (per stellar baryon) both 
decrease significantly as T e s drops below 1.5 x 10 4 K. 



2 We assume a flat spectrum across the narrow LW bands; the 
"sawtooth modulation" by the IGM has been absorbed into J cr it • 
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Figure 2. Top panel: Critical flux normalized (as in Fig. [TJ at 
the Lyman limit, as a function of effective temperature of the 
incident blackbody radiation. The corresponding masses of zero- 
metallicity stars, shown along the upp er horizontal axis, a re ob- 
tained from the data tables provided bv lMarigo et al.l j200lh . Bot- 
tom panel: The fraction of baryons that must be incorporated into 
stars in order to produce the critical flux in the top panel. 



3.3 Comparison to previous studies 

To our knowledge, the dominant importance of the radi- 
ation from low-mass stars for global star-formation was 
n ot previously discussed or qu a ntifie d, with the exception 
of Chuzhov. Kuhlen fc Shapiro! (|2007h . These authors eval- 
uated the global impact of H~ -detachment by an IR back- 
ground at later epochs, produced by recombination radia- 
tion as the IGM is becoming significantly ionized (| Glover! 
I2007I considered a similar scenario, with radiation from gas 
that is highly ionized by local sources). They showed that 
if the ionization is produced by stars with a soft spectrum, 
as we consider here, the H~ -detachment rate can be fur- 
ther boosted and can become globally important. However, 
they do not explicitly compare H~ and LW feedback, and do 
not answer the two questions addressed in this Letter, how 
many low-mass stars need to form (1) before H~ photo- 
detachment becomes the dominant feedback mode, and (2) 
for this mode to become globally important. 

Other previous works have touched on different 
aspects of the globa l radiative feedback discuss ed here . 

Mesinger. Bryan fc Haiman! (2006) 



include 



Exam ples _ . 

and iMesinger. Bryan fc Haiman! (120091 ), studying ra- 
diative feedback in a statistical sample of several 
hundred early minihalos in cosmological simulations; 



lHaimanT Abel fc Reesl <|2000T > and I Johnson? Greif fc Bromml 
12008), studying self-regulation of star-formation in mini- 
halos tlrroriglTjjadia^^ models; 
and IWhalen. Hucckstacdt & McConkid (|201Ch . studying 
radiative feedback in detailed hydrodynamical simulations 
with the stellar masses of the sources extending down 
25 M@. However, in all of these works, feedback was due to 
Ly man- Wern er (or ioniz ing U V) radiation. 



Omukail (12001 



IShang. Bryan fc Haima: 



■Bromm fc Loebl J2003h 
3 (|201Ch , IWolcott-Green fc Haiman! 



© 2012 RAS, MNRAS 000,(21] 



Infrared Radiative Feedback 5 



|201ll ') and IWolcott-Green, Haiman fe Brvanl |201ll ') have 
performed calculations similar to the one proposed here, 
including H~ photo-detachment, but have focused on 
the more massive (T v i r > 10 4 K) atomic-cooling halos. 
Because the gas in these halos can cool via neutral H and 
reach high densities, a much larger flux J cr i t is required to 
suppress H2-cooling. Nevertheless, these works have found 
results similar to the ones here: the critical flux for a soft 
spectrum is ~ (1 — 2) orders of magnitude lower than for 
a hard spectrum. F or example, using 3D simulat ions of 
three different halos, I S hang. Bryan fe Haiman] (|2010h found 
30 < J cr it < 300 (T eff = 10 4 K) versus 10 4 < J crit < 10 5 
(T e ff = 10 5 K). As found here, H~-detachment dominates in 
the former case, whereas H2-dissociation dominates in the 
latter. When self-shielding in the LW lines of H2 is modeled 
more accurately, however, this difference is reduced by about 
an o rder of magnitude (|Wolcott-Green. Haiman fc Brvanl 



4 CONCLUSIONS 

The main conclusion of this Letter is that if the mass- 
fraction of low-mass (few M Q ) stars exceeded ~ 90%, then 
their early IR background radiation dominated over the 
LW background in suppressing H2 formation. This low-mass 
fraction is comparable to those in present-day IMFs, and 
it is interesting to note that in this limit, star formation 
in minihalos could be more efficient than if the early stars 
were massive, owing to the lack of UV feedback and heating 
inside halos. The early IR background from the low-mass 
stars would then exert significant net feedback, and regu- 
late the star-formation history in the early Universe once a 
fraction /* = a few x/,,lw of baryons were converted into 
stars. The threshold /*,lw is the fraction required for strong 
LW feedback (from massive stars), which is ~ 0.3% that re- 
quired for reionization assuming n 7 = 10 ionizing photons 
per hydrogen atom. Future investigations of radiative feed- 
back in the early Universe, which include low-mass stars, 
should therefore include H~ photo-detachment. Our results 
also highlight the need for an accurate calculation of the IR 
photon output of low-mass stars. 
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